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Asteroid rotation excitation by subcatastrophic impacts 
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ABSTRACT 

Photometric observations of asteroids show that some of them are in non-principal 
axis rotation state (free precession), called tumbling. Collisions between asteroids have 
been proposed as a possible asteroid rotation excitation mechanism. We simulated sub- 
catastrophic collisions between asteroids of various physical and material parameters 
to find out whether they could be responsible for the excited rotation. For every sim- 
ulated target body after the collision, we computed its rotational lightcurve and we 
found that tumbling was photometrically detectable for the rotational axis misalign- 
ment angle (3 greater than about 15°. We found that subcatastrophic collisions are 
a plausible cause of non-principal axis rotation for small slowly rotating asteroids. 
The determining parameter is the ratio of the projectile orbital angular momentum 
to the target rotational angular momentum and we derived an approximate relation 
between this ratio and the angle (3. We also compared the limiting energy for the onset 
of tumbling with the shattering energy. Slowly rotating asteroids of diameter 100 m 
and larger can be rotationally excited by collisions with energies below the shattering 
limit. 

Key words: minor planets, asteroids: general - methods: numerical. 



1 INTRODUCTION 

Most asteroids appear to be in a basic state of rotation, 
i.e., they rotate around the shortest principal axis (the axis 
with the maximum moment of inertia). We know, however, 
that asteroids with small sizes and slow rotations are in non- 
principal axis rotation states, for which Harris ( 1994) coined 
the term tumbler. They show lightcurves with two frequen- 
cies (and their linear combinations), reflecting their rota- 
tional state of free precession. The first confirmed tumble r 
was 4179 Toutatis (jSpencer et al.lll995l . R)stro et al.lll995h . 
a small Apollo asteroid that was thoroughly observed in De- 
cember 1992 around its close approach to Earth. 

Since the non-principal axis rotation state has a 
higher energy than the basic one, we are interested in 
what processes excite asteroid rotations. Propos ed expla- 
natio ns were a solar radiation pressure toque (|Paddackl 
1969), a torque given by absorption of sun radiation and 
its re-emission as thermal radiation by irregular asteroid 
surface with a thermal inertia (the Ya rkovsky-0 'Keefe- 
Radzievskii-Paddack, or YORP, effect; iRubincaml l200d . 
IVokrouhlickv et al.l I2007T ). ejec ta escaping from t he tar- 
get asteroid after the collision (jScheeres et al.ll2000h , grav- 
itational torques during close plan etary ffybys (only pos- 
sible for planet-crossing asteroids, iRichardson et al 
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iBlack et al.l 1 1999)) . and a torque cau sed by mass ejection 
from a comet ( Peale fc L issaucr 1989). More discussion on 
this can be found in lPaolicchi et al.l ( 20021 ) . 

IPravec et al.l (|2005l ) presented a review of known tum- 
blers, their dynamical characteristics and also speculated on 
the origin of their spins. They discussed the YORP effect as 
a process that could be responsible for slowing down the 
rotation of asteroid and possibly even for their rotational 
excitation. The other possibility mentioned was a subcatas- 
trophic collision which could be much more effective for an 
asteroid previously slowed down by YORP. 

Collisions between asteroids are a fundamental 
process, affecting ast e roid 



s p ins and size distr ibution 
i|Harris fc Burns! 1 19791 . lHarrisf Il979l ; ICellino et a! 
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forming asteroid families (Hi ravamalll918l . IHiravamalll927f ) 



1991 
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and some binaries (Noll 2006), a nd also triggering activity of 
main belt asteroids ilJewitt et alj201ll . iBodewits et al.ll201ll ; 
IStevenson et al.l 120121 ). It is interesting to ask, if collisions 
can be plausible explanation for the excitation of asteroid 
rotation. 

The most common surface features of the minor bod- 
ies with rocky surface are impact craters. Some craters are 
of substantial size. Actually, the largest craters observed 
on asteroids visited by space probes have D/Rm ~ 1 and 
some even D/Rm ~ 1.5, where D is the crater diameter 
and R m is the impacted bod y mean radius — the radius of 
a sphere of the same volume i| Thomas! 1999). For instance, 
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the tumbling asteroid 253 Mathilde with the mean radius of 
26.4 km (Tcdesc o et al.||l992ri an d very slow rotation (main 
period 418 h. lMottola et al.ll 199a ), is covered by several huge 
impact craters. The largest one has a diameter of 33 km 
(|Veverkaet al.ll 19971 ). 

The process of collisional excitation of asteroid 
rotation has been exam ined in some earlier papers 
JGauchez fe Souchavll2006r ). but they did not specifically fo- 
cus on the onset of tumbling. We use a similar parametriza- 
tion of the problem, but we only assume explosive cratering 
events and concentrate on non-principal axis rotation. 

In section [5] we describe our numerical model that we 
used to study the effects of collisions between asteroids on 
their rotations, section [3] shows the dependence of a rota- 
tional axis misalignment on physical and dynamical param- 
eters of the colliding bodies, in section 2] we describe how 
the misalignment depends on the ratio of projectile's orbital 
AM to the target's rotational AM and in section [5] we plot 
specific impact energy of the collisions as a function of the 
mean target radius. 



2 SUBCATASTROPHIC IMPACT MODEL 

We investigate the effect of subcatastrophic impact on the 
rotation of an asteroid in a following way: there are two 
bodies described by a set of physical, dynamical and shape 
parameters (size, density, velocity etc.) that experience hy- 
pervelocity collision. The projectile is completely destroyed 
and an impact crater is created on the target. We then recal- 
culate the dynamical characteristics of the target, calculate 
the rotational lightcurve and check if the tumbling of the 
target can be detected by photometric observations. We de- 
scribe our model in more detail in following. 

The target is a homogeneous triaxial ellipsoid in re- 
laxed rotation state, i.e., rotating around the shortest axis. 
It is described by following parameters: the mean radius of 
the body, the semi-axes ratios a/c and b/c (a ^ b ^ c), 
the body's bulk density p, the material strength_|, initial ro- 
tation period and whether the body is porous or not (see 
chapter 12.41 for details on this). The projectile is a homoge- 
neous sphere with diameter d and the mean density 8. We 
do not take into account the rotation of the projectile since 
its angular momentum is negligible. 

We observe the impact in an inertial frame connected 
with the centre of mass of the two bodies. The frame coor- 
dinate axes are identical to the principal axes of the target 
at the moment of impact (x in the direction of the target's 
longest axis, z in the direction of the shortest one). 

The projectile is moving towards the target at a relative 
velocity v and hits it at the impact point I with spherical 
coordinates <f> (longitude) and 9 (latitude) on the surface of 
the target body. The prime meridian is defined as y = and 



x ^ 0; together with the longest target axis it makes the 
common plane. 

The impact speed is much greater than the target's 
escape velocity so we can neglect mutual gravity of the 
bodies and any curv ature of the projectile's trajectory 
(|Love fe Ahrensll 19961 ). 

The projectile is completely destroyed in the impact and 
a crater is formed on the surface of the target at the impact 
point. Its diameter and depth are calculated by using scal- 
ing laws based on point source approximatio n of the impact 
|Holsapplelll99"3l : iHolsapple fe Housed 120071) . 

The geometrical representation of the crater is a 
paraboloid of revolution; this corresponds to the observa- 
tions of simple impact c raters on the Moon and other bod- 
ies in the Solar System (|Chappelow fe Sharpton||2002| ). The 
axis of the paraboloid is perpendicular to the local sur- 
face. We do not take into account border rim of the crater 
(which is usually present in real craters); we only model 
simple craters in our simulations. Complex impact struc- 
tures can only be observed on the largest asteroids in the 
main belt and on large icy bodies in the outer Solar System 
where the gravity plays a substan tial role in crater formation 
JLeliwa-Kopystvriski et aLll2008l ). 

The diameter of the largest crater that a body of a 
gi ven size can bear is calculated accord ing to the relation 
in lBurchell fe Leliwa-Kopystvnskil (|2010l ). They studied the 
maximum crater size on small bodies (asteroids, icy satel- 
lites, comets), which were imaged by space probes, and 
showed, that it follows the relation 



D = -(0.17 ± 0.10) + (1.01 ± 0.08)iC 
0.7 km < Rm < 120 km, 



(1) 



D being crater diameter and i? m the target mean radius 
both in kilometres. 

The relative diameters of the impact craters F — 
D/Rm, formed on the target in our simulations, were usually 
kept less than 1.2 6, which is the value for the larg est crater 
on 253 Mathilde (|Leliwa-Kopvstvriski et al.l 120081 ). In some 
simulations we allowed for higher values of up to F — 1.6, 
since there is indirect e vidence of large crater on 90 Antiope 
IjDescamps et al.ll2009l ). However, conclusions based on this 
higher limit should be made with caution. 

For every collision we calculated also its specific im- 
pact energy (kinetic energy of the collision per unit mass 
of the target ) and compared it to t he di spersal criterion as 
described bv lStewart fe Leinhardtl (|2012f) . We plot our sim- 
ulation results in a graph of specific impact energy vs. mean 
target radius and discuss the results in section [5] 

The ejecta are assumed to leave the impact point with 
a cone-shaped velocity field symmetric around the line per- 
pendicular to the surface at the impact point. We take this 
effect into account when calculating the angular momentum 
transfer efficiency. 



1 In most of our simulations, material constants used involve 
some microporosity. In some simulations we modelled macropo- 
rosity as well. 

2 We assume effective strength for the excavation, because the 
shock wave fractures the body during the early phase of the im- 
pact (it is not tensile, shear or other materi al strength that ca n 
be measured in la bora tory); for details see iRichardsonl (2009), 

iNolan et all J1996I ) and iHolsappiel J2007I) . 



2.1 Inertia tensor calculation 

The inertia tensor of the target before the impact is calcu- 
lated analytically for the body being homogeneous triaxial 
ellipsoid. The inertia tensor after the impact is calculated 
numerically since the body is irregularit is the triaxial ellip- 
soid with the crater on its surface. 
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The body surface is divided into triangular facets and 
together with the origin they form the tetrahe dron. Its iner - 
tia tensor is calculated analytically (see, e.g.. iTononl 12004 ) 
and using the Steiner's parallel axis theorem it is translated 
into the original centre of mass coordinate system (reference 
frame, origin O). The inertia tensors of all such tetrahedra 
are then summed up and the inertia tensor of whole body is 
obtained. Then it is translated to the new reference frame 
with origin O* , since the centre of mass has shifted due to 
the impact. 

The situation with facets covering the crater is only a 
bit more complicated. First the crater border is calculated. 
Then all the facets' vertices of the ellipsoid surface inside this 
circular border are replaced with points of the paraboloid 
describing the crater floor. New facets' vertices are made of 
these points and the same procedure of tetrahedra inertia 
tensor calculation is used then. 



2.2 Principal axis deviation calculation 

We try to find out whether a subcatastrophic impact can 
create an observable excitation of rotation of the impacted 
body. The target's rotation was relaxed before the impact. 
It means it was rotating about the shortest axis, i.e., the 
axis with the largest moment of inertia for the homogeneous 
ellipsoid. We denote this moment of inertia 73 and it holds 
that h < h < Li- 

An indication of the excited rotation will be a nonzero 
angle between the rotational angular momentum vector of 
the body after the impact L* and its shortest principal axis 
with the largest moment of inertia 1$. The procedure of its 
calculation is as follows. 

First the inertia tensor of the target after the impact 
is calculated. Since the body possesses the impact crater, it 
lacks symmetry and the inertia tensor is calculated numeri- 
cally. 

Because the centre of mass of the target O* is shifted 
with respect to the original centre of mass O, the resulting 
inertia tensor has to be shifted accordingly (the method cal- 
culates it with respect to O, hence it is not a central inertia 
tensor for the new body). We want to have the inertia tensor 
in a body frame (which is not inertial) with an origin in its 
centre of mass, because then the components of this inertia 
tensor do not change with time as the body rotates. The 
transformation reads 



vector just after the impact. We use the momentum conser- 
vation law 



l*fc = Uk + M(a Sik — a,ia k ) 



(2) 



where M is the total mass of the body, a is the centre of 
mass shift vector (between O and O*), a is its magnitude, 
ai and at are its cartesian components and Sik is Kronecker 
delta (see, e.g.. lKvas nica 1997). 

Now we have to find the principal axes of the target after 
the impact. We can calculate them by finding the diagonal 
form of its inertia tensor, because they are determined by 
the eigenvectors of this tensor. The diagonalization is given 
by 

Id = Q T I*Q, (3) 

where Id is diagonal inertia tensor, Q and Q T are orthogonal 

matrix and transposed orthogonal matrix, respectively. Its 

rows (and columns, respectively) are the eigenvectors of I*. 

In the next step we calculate the angular momentum 



m t vt + m p v imp — m v 



eject a •, 



and angular momentum conservation law 

L>t + -L/p + -L/orb,t + -Lorb,p = J-/ + -L or b H~ -dejecta , 



(4) 



(5) 



where m p is the mass of the projectile, «i mp is the projectile 
velocity vector, which for very low mass ratio of the projec- 
tile to target approaches the impact velocity, m t and m* 
are the masses, vt and v* are the velocities of the target be- 
fore and after the collision, respectively, and P e jecta is the 
momentum of the ejecta. Lt and L p are the rotational an- 
gular momenta of the target and the projectile before the 
impact, respectively, L rb,t and L rb, P are the orbital an- 
gular momentum vectors of the target and the projectile, 
respectively, with respect to the system's centre of mass, L* 
and L* rb are the rotational and orbital angular momenta 
of the target after the impact, respectively, and Lejecta is 
the angular momentum of the material ejected during the 
impact. Since the projectile has dissintegrated during the 
impact, all the quantities on the right sides of Eqn.|3]and[5] 
describe the postimpact target body, and so we omit the 
subscript t . 

The orbital angular momenta are calculated as 



^orb,p — Hlp\Xp X V\ m p) 

iorb.t = m t (r t X v t ) 
A-orb = m (r xv) 



(6) 



where r p , rt and r* are the radius vectors of the projectile 
and the target before and after the impact, respectively, with 
respect to the system's centre of mass. We calculate the 
target velocity before the impact as 



v t 



nit 



(7) 



Without an appropriate model of the ejecta, we need to 
at least roughly estimate its dynamical effect. We decided 
to use the efficiency of the linear (and angular) momen- 
tu m transfer measur e d in l abora tory impact experiments 
by lYanagisawa et al.l (|l996T ) and lYanagisawa fc Hasegawal 
(2000) and we assumed that this relationship holds also for 
large scale impacts. Using the efficiencies rj and f, the lin- 
ear momentum transfer efficiency in normal and tangential 
direction with respect to the local surface, respectively, we 
calculated the v* including the effect of ejecta. For the def- 
inition of these efficiencies and more detailed discussion see 
section 12.31 

Lastly we needed to describe Lejecta. For small water- 
ing impact we can calculate the difference L or b,p — -dejecta 
as 

L QT b,p — iejecta = 1)171? (Vp X Unorm) + C m P (^P X 1>tang) , (8) 

where u norm and w t an g are impact velocity components nor- 
mal to and tangential to the local surface of the target, re- 
spectively. 

The rotational angular momenta of the bodies are sim- 
ply given by 

it = lt«t and Lp — \ p u} p , (9) 

where we introduced inertia tensors of the target l t and of 
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the projectile l p and angular velocities ojt and u> p of the 
target and the projectile, respectively. 

When we put the equations[Bl[71[S]and[5]to the equation 
of angular momentum conservation O we can calculate the 
angular momentum vector of the target after the impact as 

L — ltU>t + IpkJp + (iorb.p — iejecta) + iorb,t _ Lorb -(10) 

Finally we can calculate the angle between the rota- 
tional angular momentum vector of the target after the im- 
pact L* and the axis with the largest moment of inertia I£ 
(the axis is represented by unit vector E3, the eigenvector 
of the inertia tensor I*). We denote this angle as /3, 



/3 — arccos 



L E3 

\L*\\E 3 \ 



(11) 



The angle /3 is not constant, but it varies with changes 
of the principal vector E3 as the body precesses. We calcu- 
late the value of /? right after the collision and we use it as 
an approximate measure of the magnitude of tumbling dis- 
played in lightcurve, for which the use of the instantaneous 
value of the misalignment angle is sufficient. 



2.3 Angular Momentum Transfer Efficiency 

During the collision the projectile angular momentum is 
transferred to the target. The target body, however, does not 
obtain all of the angular momenta carried by the projectile, 
as shown by the experime n ts of small scale impacts (s e e, e.g. , 
lYanagisawa et al.l (|199a ). lYanagisawa fc Hasegawal (2000) 
and references therin). A part of the angular momentum is 
carried away by high velocity ejecta and spall fragments and 
so we introduce the momentum transfer efficiencjo defined 



V = 



Wlp 'Wnorm 



and 



c = 



771 ^tang 
Wlpvtatig 



(12) 



where v* is the target speed after the impact, v is the projec- 
tile preimpact speed. The subscripts norm and tang denote 
the normal and tangential components to the target's local 
surface at the impact point, respectively. In our calculation, 
we use the experimentally obtained values of the efficiencies, 
but the thing is actually a bit more tricky. 

From the impact experiments we see that the efficien- 
cies depend on impact angle a, which is an angle between 
the impact velocity vector and the line perpendicular to the 
surface at the impact point. The experimental results give 
the following dependencies 

n(a) = 1 + 770 cos a and £(cr) = fo cos a . (13) 

lYanagisawa fc Hasegawal (|2000P ) obtained the values of con- 
stants 770 and £0 at 4kms~ 1 impact speed, which is close 
to the media n impact speed bet ween the main belt bodies 
(~ Skms -1 . iBottke et al.l I2005T ) . For basaltic targets they 
give 770 = 1.52 and (0 = 0.409. For sand target which could 
be similar to a regolith-covered asteroid. lYanagisawal (2002) 
gives £0 = 0.687. He did not obtain the 770 value and there- 
fore he uses 770 = (77 = 1) for an asteroid large enough to 

3 The efficiency can be larger than one because backfire ejecta 
carries away some momentum in the direction opposite to the 
projectile course. 



reaccumulate all the ejecta by its gravity and we follow this 
approach. 

There is also a relation between th e efficiencies a nd th e 
shape of the target body as described in lYanagisawal |2002). 
It is sufficient to consider the equations given above since 
the impact angle a is different for a sphere and an ellipsoid 
for the same impact velocity and therefore we get different 
values of 77 and C, for different target shapes. 



2.4 Scaling laws 

To find out how large will be the cra ter formed on th e as- 
teroid, we followed the scaling laws of Holsapplc (2003) and 
iHolsappld l| 19931 ). The scaling laws in general predict the 
outcome of an experiment from another one with different 
values of problem parameters. When we study impact pro- 
cesses, we have centimetre-scale impact experinrents from 
the lab and we try to predict the outcome of the same ex- 
periment on a much larger (kilometre) scale. 

Based on shapes of simple craters on the Moon and 
other Solar System bodies we decided to model the impact 
crater as paraboloid of revolution. We assume that the af- 
termath of the crater formation (e.g., a collapse of the crater 
walls) does not have a substantial dynamical effect (the mass 
redistribution is small). 

The radius and depth of the crater is given by the fol- 
lowing formulae 



Re 



K X V, 



1/3 



and 



D c 



K d V, 



1/3 



(14) 



where K r and Kd are the shape con stants specific for a ma- 
terial of the asteroid (|Holsapplell2003i ) and V c is the volume 
of the crater given by 



4nd 3 5 

Vc = — 7T V . 

3 p 



(15) 



d being the impactor diameter, 5 its density, p is the target's 
density and 7ry is given by 



ttv = K 1 I 1Y 2 



(I) 



(6i>-2-n)/3/i 



K 2 



7T:i 



(?) 



(6f-2)/3fi 



(2+M)/2 



-3m/(2+m) 



(16) 



7T2 



gd 

U 2 ' 



7T 3 = 



Y 

~plP 



Constants K\ and K 2 and expo nents u and v c ome from 
laboratory impact experiments JHolsapplel ll993T ). g is the 
gravitational acceleration at the impact point, U is the im- 
pact velocity component perpendicular to local surface and 
Y is the target's material strength. In our simulations we 
consider two cases, a rocky material corresponding to non- 
porous bodies and lunar regolith for porous asteroids. For 
the values we use, see Table [I] 



2.5 Synthetic lightcurves 

We generated rotational lightcurves for the resulted nrodel 
bodies. We used th e program written by Mikko Kaasalainen 
|Kaasalainenl I2001T ) and modified by Josef Durech (2011, 
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Table 1. Scaling and material constants for non-porous and porous body materials. 



Material 


K x 


K 2 


M 


V 


Y (MPa) 


p(gcm 3 ) 


K r 


K d 


Rock 
Lunar rcgolith 


0.095 
0.132 


0.257 
0.26 


0.55 
0.41 


0.33 
0.33 


10 
0.01 


2.0 
1.5 


1.1 
1.4 


0.6 
0.35 



Table 2. Physical parameters of the sample target bodies. 



principal axis rotation with minor deviations 



Figure 



target mean 
diameter [km] 



initial rotation 
period [h] 



[dcg] 



1 


144.2 


32 


1.2 


m 


36.1 


32 


11.5 


3 


1.4 


2 


15.5 


HI 


7.2 


12 


18.8 


m 


7.2 


22 


34.3 



principal axis rotation 




0.4 0.5 0.6 
rotational phase 

Figure 1. A lightcurve of a body impacted by a small projectile, 
/3 = 4.2°. The horizontal axis is the rotational phase, the vertical 
axis is a light flux normalised to unity. It is undistinguishable 
from a lightcurve of principal axis rotator. 



pers. communication), which generates the lightcurve for an 
arbitrarily shaped body with given dynamical parameters. 
The shape is given by the set of oriented triangles building 
up the bo dy surface mesh, th e light scattering law was that 
of Hapke (jBowell et al.l ll989h The input also includes the 
orientation of axis of rotation at some instant of time as 
well as the Euler angles for that instant, that determine the 
body's orientation in space. 

The Euler equations of motion are then solved and the 
integral light flux is calculated for given observational ge- 
ometry with a chosen time increment. We simulated real 
observations with varying time sampling rates to avoid alias 
effects in generated lightcurves and we also added a Gaus- 
sian noise. 

A blind test lightcurve analysis was done for a set of 
cases with increasing /? to find out if the tumbling is de- 
tectable in the synthetic lightcurve. In Figs.[T}{5]we present 
five sample lightcurves. The axial ratios of the target were 
a/c — 2.0 and b/c = 1.5. The diameter of a projectile was 
about 5% of the mean diameter of the target body in all 
these simulation runs. Other parameters describing the body 
are in Tabled 
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Figure 2. A lightcurve with apparent minor deviations from sin- 
gle periodicity, /3 = 11.5° , only marginally distinguishable by pho- 
tometric observation. 
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0.9 1 



Figure 3. A moderate tumbling signal is apparent in this 
lightcurve, [3 = 15.5°. 



The values of rotation axis misalignment for which the 
tumbling is clearly detected start at ft ~ 15° . In section [4] 
we will show that the value of /3 is closely related to L or b/Lt 
ratio independent of other parameters describing the two 
bodies. 



3 ROTATION AXIS MISALIGNMENT 

In this section we present the results of our simulations. 
In individual runs, we varied one or two input parameters 
(e.g., body size or its rotation period) while keeping other 
constant. Then we plotted graphs showing the relationship 
between the varied parameters and the misalignment /? of 
the angular momentum vector after the collision (see section 



© 2013 RAS, MNRAS OOO,[Q{T0] 



6 Tomds Henych, Petr Pravec 



tumbling 



120 



1.3 



1.1 



I I I I I I 




•****». 


""», 


/v ••••. 


.- -". 


« * * 






• v ♦ •» Sk * • i ♦ 










- 0.9 - 

0.8 

0.7 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

rotational phase 
Figure 4. In this lightcurve the tumbling is evident, /3 = 18.8°. 
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Figure 5. Obvious tumbling for this slowly rotating body, 
34.3°. 



I2.2I) . The following graphs do not show results of all our 
simulations, but rather a representative sample. 

The constant parameters describing the body in these 
runs were a/c = 2.0, b/c = 1.5, Y = lOMPa, bulk densities 
of target and the projectile were 2 gem -3 . The projectile 
collided with the body at v x = — 5kms _1 and the impact 
point coordinates were (j> = 15° and 9 = 15°. The incidence 
angle for such geometry is approximately 35°, which is ten 
degree s less than the statisti cal average for a spherical target 
body (|Love fc Ahrenslll996ft . 

Figure [6] shows the misalignment of the AM vector for 
various values of the initial rotation period and four body 
sizes. As expected, /3 increases with increasing rotation pe- 
riod. The slower rotation of the target implies it has lower 
angular momentum L t and therefore the projectile with a 
given orbital angular momentum L or b will excite the body 
more. 

We can clearly see, that even for the largest body the 
tumbling could be detected for periods longer than about 
42 h. To understand qualitatively this graph, we can use di- 
mensional analysis. For the rotational angular momentum 
of the target and for the orbital angular momentum of the 
projectile, respectively, we have 
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Figure 6. A relationship between the AM vector misalignment 
and the initial rotation period of the body, for four different sizes. 
The projectile size was adjusted so that it excavated a crater of 
maximum size (see section [2j. 



T T 2 5 27T 



(17) 



where rt is the mean target radius and Pt is the target's 
period of rotation. The ratio of these two is 



Lorb 



,.4 



(18) 



It means that the smaller the target and the longer the pe- 
riod of rotation, the higher the ratio of the angular momenta 
and therefore higher /3. 

In another run, the sampled parameter was the projec- 
tile size. The results are plotted for small body both porous 
and non-porous (Fig. [7|. As the relative projectile size in- 
creases, /3 also increases. A larger projectile makes a larger 
crater on the target body and hence larger change of the 
inertia tensor. 

The parameters describing the body and the impact 
were the same as previously, the mean body diameter was 
1.4 km. For a porous target, there is a different set of the scal- 
ing and crater shape parameters as described in section [2~4l 
and different cratering efficiency (see section |2~3]) . There are 
only subtle differences of at most two degrees in /3 between 
porous and non-porous targets. 

In Fig. [HI we plot the strength-misalignment relation- 
ship for a small and a moderate size body. For increasing 
target strength there is an increasing threshold crater size 
(see section[2]for details). We kept the actual crater size near 
the threshold, so for larger target strength we had larger 
projectile with larger L or b (Lt being constant). 

Other parameters than the strength and the body sizes 
(144 m and 7.2 km, respectively) were the same as in previ- 
ous experiment, the initial rotation period was 32 h. 

The last run shows a dependence of the misalignment 
on the target size, see Fig. [9] It is clearly seen that with 
decreasing target size /3 increases. As before, it was caused 
by increasing L or b/L t ratio, which can be, again, explained 
by dimensional analysis (see Eqn. ll8[) . 
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Figure 7. A relationship between the AM vector misalignment 
and the projectile size for small target. The simulation for non- 
porous and porous body. 
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Figure 8. A relationship between the AM vector misalignment 
and strength of the small and moderate size target. 
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Figure 9. A relationship between the AM vector misalignment 
and target size. 
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Figure 10. All the previous results included in one graph, /3 as 
a function of L or b/£t- 



4 ANGULAR MOMENTA RATIO 

Here we show, that the results of all our simulations can 
be described by a ratio of the projectile orbital to target 
rotational angular momentum L OY ^/L t . The relationship can 
be approximately described by the formula 



cos/3 = ± 



sin ip 



(L t /L olh -fcosVO 2 



-1/2 



(19) 



where the + sign is for L t ^ — L or b cost/; and — sign is for 
Lt < — iorb cos ip, ip is the angle between the two angular 
mometum vectors just before the collision. For derivation of 
this formula see Appendix lAl 

We plotted all the simulation results in Fig. 1 101 together 
with the curve given by Eqn. 1191 

There is another obvious feature in the above mentioned 
graphthe values of (3 converge to the value of ip. This is to 
be expected when the ratio of the magnitudes of these two 
vectors becomes much greater than one. 

We have checked this convergence also for other an- 
gles (other impact geometries), the results are presented in 
Fig. 1111 The target was 50 m body with other parameters as 
in previous simulations and the projectile collided with it at 
V x = — 5kms _1 . The impact point coordinates are described 
in Fig. [TT]by (j> and 9 (see section 0). 



5 DISPERSAL AND SHATTERING CRITERIA 

It is important to compare the specific energy of the col- 
lision with the shattering or dispersal criterion so that the 
collision is subcatastrophic (cratering regime). Therefore we 
plotted specific impact energy as a function of the mean 
target radius and calculate the shattering criterion for most 
collisions we simulated. We considered targets with various 
dimensions, but constant a/c = 2.0, b/c = 1.5, with a ro- 
tation period of 36 h, mean density 2 gem -3 and strength 
1 MPa. The projectiles were of various sizes, their density 
was 2 gem -3 and the impact velocity was v x — — 5kms _1 . 
The impact point coordinates were <p = 21° and 9 = 21°, 
the incidence angle for such geometry is approximately 44°. 

The results are shown in Fig. 1121 

In the calculation, we followed IStewart fc Leinhardtl 
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Figure 11. Test of convergence to ij> for various impact geome- 
tries, or angle between Lt and L or b vectors. 



(|2012f ) who derived the dispersal criterion in gravity regime 
from their set of numerical simulations of collisions for wide 
variety of input parameters describing the colliding bodies. 
As shatterin g criterion we take 1/4 of the dispersal value 
according to Irlousenl (J2009J)- Since the diversity of disper- 
sal criteria in literature is rather high, we plot other such 
functions in Fig. [T^] for comparison to our choice. 

We also plot the limiting value of the specific energy for 
which the tumbling after the collision becomes detectable in 
the lightcurve as described in section 1231 For its derivation 
see Appendix iBl The bodies, for which the specific impact 
energy of the collision was less than this limiting value, are 
labeled 'not tumbling' for simplicity, but this label implies 
that the tumbling would not be detected in their lightcurves. 
If the collision energy exceeds the shattering criterion, we 
label the resulting body 'tumbling but shattered' and we 
mark it with different symbol (x) in the graph. 

We plotted an upper limit of the specific impact energy 
in the graph. It corresponds to the ratio of the diameter of 
the largest crater to the mean target radius D/R m — 1.26 
(the largest known crater on 253 Mathilde). 



6 SUMMARY 

We studied effects of subcatastrophic collisions on asteroid 
rotation. We found that it is plausible that such collisions 
are responsible for observed excited rotation of the small 
bodies. The important parameter of the collision is the ratio 
of the projectile orbital angular momentum to the target 
rotational angular momentum before the collision (L OT b/Lt). 
We found an approximate relation between the rotational 
axis misalignment angle and this ratio, which can be used 
for estimation of the excitation effect of a particular collision 
between two bodies. 

We plotted the results of our simulations to the graph of 
the specific impact energy vs. mean target radius (Fig. [12]). 
We showed that for a reasonable choice of input parameters 
and the dispersal limit, the excitation of rotation can take 
place below the shattering limit for a target body. 

In a future work, we plan to estimate the probability of 
subcatastrophic collision capable of rotational excitation for 
bodies of various sizes during their lifetime. Also we will ac- 
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Figure 12. The specific energy of collision as a function of the 
mean target radius. Small red circles denote bodies for which 
the tumbling was not apparent in the lightcurve. Large green 
circles are tumbling bodies which did not exceed the shattering 
criterion. Green cross symbols are for bodies that were shattered 
or even had a half of its mass di spersed by the collision. Th e 
dispersal criterion is according to IStewart &: Leinharda (2012). 
Other dispersal criteria are also displayed. The thick black line is 
the limiting specific impact energy as derived in Appendix 151 



count for the damping timescale of tumbling asteroids to find 
out the probability of finding such bodies in non-principal 
axis rotation state. 
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APPENDIX A: MISALIGNMENT 
ANGLEANGULAR MOMENTA RATIO 
DEPENDENCE 

Equation [11] defines the rotational axis misalignment angle 
after the collision. As we showed, it is closely related to the 
ratio of L or b/L t . Here we derive the approximate relation 
between these two quantities. 

We assume that the shortest principal axis (represented 
by vector E3) of the target remains unchanged during the 
collision (the small deviation is caused by the crater formed 
on the body). We only consider target rotational angular 
momentum Lt and projectile orbital angular momentum 
£orb, P (further only L or b for simplicity) and neglect all 
other angular momenta in Eqn. [5l so the resulting L* is 



L — -Lorb + Lt 



cos/3 = 



Hence we have 

L*E 3 _ (L rb + Lt) 



\L*\\E S \ 



I Lorb + Lt\ 



(Al) 



(A2) 



where the subscript z denotes the z component of the vector. 
This can be further written as 



cos/3 



(LorbCOS?/' + L t ) 



[(L orb cos V + L t ) 2 + L 2 ovh sin 2 V] 1 / 2 



(A3) 



where ip is the angle between _L or b and Lt vectors. After 
some algebra we have 



cos/3 = ± 



1 + 



sin 2 ip 



(Lt/Lorb +cosV0 2 



-1/2 



(A4) 



where the + sign is for L t ^ —Lorb cosip and — sign is for 
Lt < — L or b cos ip. 



APPENDIX B: LIMITING SPECIFIC IMPACT 
ENERGY 

Here we derive the limiting value of the specific impact en- 
ergy Q tumb as a function of physical parameters of the col- 
liding bodies and the ratio of their angular momenta (orbital 
angular momentum of the projectile L or b and the rotational 
angular momentum of the target before the collision L t ). 
The ratio of the angular momenta is given by 



L t _ 2n(a 2 + b 2 )abcp 
L rb 5r^Vim P SPtr t 



(Bl) 



where Ui mp is the impact speed, 5 and p are the projectile 
and target density, respectively, n and r p are their mean 
radii, respectively, a, b and c are the target's semi-axes and 
P t is the target's period of rotation. The specific energy of 
the collision in the system's centre of mass for m p -C nit is 



m P^imp ^p^imp^ 



2m t 



abcp 



(B2) 



mt and m p being the target and the projectile mass, respec- 
tively. 

We solve Eqn. lA4l for L t /L or b ratio, Eqn. lBll for r p and 
insert both of them into Eqn, IB2l We then obtain the specific 
energy for which the tumbling starts to be apparent in the 
lightcurve 
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n tumb 7t(a 2 + b 2 )v imp ( sin 2 ip 2 



cos^ tp . (B3) 



5rtPt V cos 2 P — 1 

In Fig. [12] we use the derived value of /? = 15° 



This paper has been typeset from a TffjX/ F/TfrjX file prepared 
by the author. 
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